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Erosion of rocky coasts spontaneously creates irregular seashores. But the geometrical irregularity, 
in turn, damps the sea- waves, decreasing the average wave amplitude. There may then exist a mutual 
self- stabilisation of the waves amplitude together with the irregular morphology of the coast. A 
simple model of such stabilisation is studied. It leads, through a complex dynamics of the earth-sea 
interface, to the appearance of a stationary fractal seacoast with dimension close to 4/3. Fractal 
geometry plays here the role of a morphological attractor directly related to percolation geometry. 
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Coastline morphology is of current interest in geophys- 
ical research [lj , and coastline erosion has important eco- 
nomic consequences Q- Also, the recent concern about 
global warming has increased the demand for a better 
understanding of coastal evolution. At the same time, 
although the geometry of seacoasts is often used as an 
introductory archetype of fractal morphology in nature 
0, Q , there has been no explanation about which phys- 
ical mechanism could justify that empirical observation. 
In the field litter atureQ one can read: "As a matter of 
some urgency, researchers concerned with coastal evolu- 
tion should consider the alternative models, even if there 
are few supporting data. The ideas of ... stochastic de- 
velopment, . . . and criticality, all deserve investigation. " 
The present work propose a minimal, but robust, model 
of evolution of rocky coasts towards fractality. 

The model describes how a stationary fractal geom- 
etry can arises spontaneously from the mutual self- 
stabilization of coast morphology and sea eroding power 
@|. If, on one hand, erosion generally increases the geo- 
metrical irregularity of the coast, on the other hand this 
increase creates a stronger damping of the sea and a con- 
sequent diminution of its eroding power. The increased 
damping argument relies on the studies of fractal acous- 
tical cavities, which have shown that viscous damping is 
augmented on a longer, irregular, surface 0,0- In the 
following, a minimal two-dimensional model of erosion is 
introduced which leads to the spontaneous evolution of 
a smooth seashore towards fractality as shown in Fig. 

Rocky coasts erosion is the product of marine and at- 
mospheric causes 9]. There exist many different erosion 
processes: wave quarrying, abrasion, wetting and drying, 
frost shattering, thermal expansion, salt water corrosion, 
carbonation, hydrolysis. A simplified picture is used here 
by assuming that the different processes can be separated 
into two categories: "rapid" mechanical erosion (namely 
wave quarrying) and "slow" chemical weakening. The 
justification is that mechanical erosion generally occurs 
rapidly, mainly during storms, after rock has been al- 
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FIG. 1: Time evolution of the coastline geometry. Upper left: 
situation at the beginning of the erosion process. The coast 
is irregular but not yet fractal. Upper right: the seashore 
at the end of the "rapid" mechanical erosion. Lower left and 
right: situations after 10 and 20 runs of the "slow" weathering 
process (with e = 0.01). 

tered and weakened by the slow weathering processes. 
We first study the supposedly rapid erosion mechanisms. 
Then we show that the full complex dynamics, involving 
fast and slow processes, build a dynamic equilibrium that 
preserves the fractal shape of the coast. 

The sea, together with the coast, is considered to con- 
stitute a resonator. It is assumed that there exists an 
average excitation power of the waves Pq. The "force" 
acting on the unitary length of the coast is measured by 
the square of the wave amplitude \I/ 2 . This wave ampli- 
tude is related to Pq by a relation of the type \I/ 2 ~ PoQ 
where Q is the morphology dependent quality factor of 
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FIG. 2: Illustration of the erosion process. The thick num- 
bers at the square centres represent the lithology {xi}. The 
numbers in the corners are the corresponding resistances n 
which depend on the local environment as explained in the 
text. The sites marked with 1 in the lower corner are earth 
sites with no contact with the sea. On the left and on the 
right respectively the situations before and after a single ero- 
sion step with f{t) = 0.5. After this step resistances are 
updated due to the new sea environment. 



the system; the smaller the quality factor, the stronger 
the damping of the sea-waves. There are several causes 
to waves damping. Since the different loss mechanisms 
occur independently, the quality factor satisfies a relation 
of the type 
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where Q coast is the quality factor due to the viscous dissi- 
pation of the fluid moving along the coast and the nearby 
islands and Qother is related to other damping mecha- 
nisms (e.g. bulk viscous damping). Studies of fractal 
acoustical cavities 0, have shown that the viscous 
damping increases roughly proportionally to the cavity 
perimeter. Therefore, one can, in first approximation, 
assume that Q CO ast is inversely proportional to the coast 
perimeter L p (t) whereas Qother is independent of the 
coast morphology. In other words, the sea exerts a ho- 
mogeneous erosion force f(t) on each coast element pro- 
portional to ^ 2 (t): 
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where L p (t) is the total length of the coast at time t (then 
L p (t = 0) = Lo). The factor g measures the relative con- 
tribution to damping of a flat shore as compared with 
the total damping. Since g represents the importance of 
coast dissipation with respect to the other mechanisms, 
which are thought to be dominant, computations are per- 
formed under the condition g <C 1. The quantity fo is 
the renormalized value of Po such that f(t)<l at all t. 

The mechano-chemical properties of the rocks consti- 
tuting the coast, which are linked to structure and com- 
position defining their "lithology" , are unknown and ex- 
hibit some dispersion. The "resisting" random earth is 




FIG. 3: Time dependence of the erosion force f(t). The 
figure shows the evolution of the sea erosion force acting on 
the coastline, for f{t = 0) = 0.65 and different values of the 
scale gradients g. The dynamics spontaneously stops at a 
value weakly dependent on g. Data refer to systems with 
Lo = 1000. Each curve is obtained by averaging over ten 
different realisations of the earth lithology distribution. 

then modelled by a square lattice of random units of 
width Lo. Each site represents a small portion of the 
earth. The sea acts on a shoreline constituted of cells, 
each one characterised by a random number x^, between 
and 1, representing its lithology. The erosion model 
should also take into account the fact that a site sur- 
rounded by the sea is relatively weakened as compared 
with a site surrounded by earth or other coast sites. 
Hence, the resistance to erosion T{ of a site depends on 
both its lithology and the number of sides exposed to the 
action of the sea. This is implemented through the fol- 
lowing weakening rule: sites surrounded by three earth 
sites has a resistance T{ = X{. If in contact with 2 sea sites 
the resistance is n = x\ . And, if site i is attacked by 3 
or 4 sides, it has zero resistance. The iterative evolution 
rule is simple: at time t, all coast sites with ri < f(t) are 
eroded, and then L p (t) and f(t) are updated together 
with the resistances of the earth sites in contact with 
the sea. Then, from one step to the next, some sites are 
eroded because they present a "weak lithology" while 
some strong sites are eroded due to their weaker stability 
due to sea neighbouring. An example of local evolution 
is shown in Fig. 

In the first steps of the dynamics, the erosion front 
keeps quite smooth and it roughens progressively as 
shown in Fig. ^ upper row. During the process, finite 
clusters are detached from the infinite earth, creating 
islands. At any time, both the islands and the coast- 
line perimeters contribute to the damping. As the total 
coastline length L p (t) increases, the sea force becomes 
weaker. At a certain time step £/, the weakest point of 
the coast is stronger than f(tf) and the dynamics stops. 
This shows that erosion reinforces the coast by preferen- 
tial elimination of its weakest elements until the coast is 
strong enough to resist further erosion. The time evolu- 
tion of f(t) is shown in Fig. 03 The fact that the process 
spontaneously leads to a total stop of erosion for a non- 
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FIG. 4: Box-counting determination of the coast fractal di- 
mension. The straight line is a power law with a slope -4/3. 
The data refer to a system with Lq — 10 4 , and g — 10 -4 . 
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FIG. 5: Scaling behaviour of the coast width a. The straight- 
line is a power law fit with the GP exponent -4/7. Each point 
is an average over 400 samples with Lq = 5000, f(0) = 0.55. 



zero wave power is remarkable. It is a direct consequence 
of random percolation [l(| as discussed below. At time 
f(tf) the coastline is fractal (see Figs. [I] and 21 ) with 
dimension Df = 1.332(3), very close to 4/3, up to a char- 
acteristic width o~. This width a is defined as the stan- 
dard deviation of the final coastline depth ((y — Y) 2 ) 
(where Y = (y) is the average erosion depth). Fig. 0] 
shows the box-counting determination of the fractal di- 
mension. The fact that Df 4/3 confirms the relation 
with percolation, as 4/3 is the dimension of the so-called 
accessible percolation cluster 11]. Note that this value 
is also close to the observed fractal dimension for many 
seacoasts, as the United States eastern shore . 

A more detailed study of the coast width a suggests 
that the model falls in the universality class of Gradient 
Percolation (GP) As shown in Fig. [51 the coast 

width o~ follows a scaling law with respect to g 

a ~ g~ a ° (3) 

with a a ~ 4/7. This law is characteristic of GP where 
sites are occupied at random with an occupation prob- 
ability that varies in space along one fixed direction be- 
tween 1 and with a constant gradient g, and a is the 
width of the frontier of the infinite cluster. The fact that 
here g is proportional to a gradient of occupation prob- 
ability is explained by the following argument. At time 
t, the erosion power is f(t) while the sea has eroded the 
earth up to an average depth Y(t), an increasing func- 
tion of t. Inverting this function, / can be written as 
f(t(Y)). There exists then a spatial gradient of the oc- 
cupation probability by the sea. For small enough g one 
can write \df/dY\ = \(g/L )(dL p (Y)/dY)\. The quan- 
tity dL p (Y) j dY is a function of g but to the lowest order 
it is a constant independent of g since even with g = 0, 
there will be an erosion due to randomness and a conse- 
quent perimeter evolution L p (t). Then to lowest order, 
the real gradient df / dY is linear in g, which can then be 
called the "scale gradient". Note that the formation of 



a fractal interface due to the spontaneous appearance of 
a gradient has already been observed in the corrosion of 
an aluminium film pjj Q3 . 

Of course, the real dynamics of the coasts are more 
complex than the "rapid" process considered above. 
They result from the interplay with the slow weather- 
ing processes, generally attributed to carbonation or hy- 
drolysis. These processes act on longer, geological, time 
scales. In order to simulate this evolution, the lithology 
parameter Xi of all the coast sites is decreased by a small 
fraction e, i.e. x\ = (1 — e)xi with e«l after the erosion 
has stopped at tf. One or a few coast sites then become 
weaker than f(tf) and the rapid erosion dynamics starts 
again up to a next arrest time. This procedure is then 
iterated. At each restart of erosion, a finite portion of the 
earth is eroded. The repeated effect of the "slow" weath- 
ering mechanisms gives a fluctuating behavior of f(t) and 
generates a coastline erosion drift with fluctuations or 
avalanches, but it does not alter the fractal features of 
the coast (see Figure 1 lower row). In that sense fractal 
geometry plays the role of a statistical attractor. In the 
language of coastal studies [13 , the system state evolves 
through a dynamical equilibrium where small perturba- 
tions (even without special trigger events) may stimulate 
large fluctuations and avalanche dynamics. This is due 
to the underlying criticality of percolation systems. In 
the authors' mind, it is this stationary regime which cor- 
responds to the geomorphologic observations of fractal 
seacoasts. 

The above simple model presents limitations, but the 
simplifications that keep the system in the universality 
class of percolation are unimportant. The use of different 
weakening rules would eventually modify the dynamics 
but not the final fractal dimension. The time separa- 
tion between rapid erosion and slow chemical weakening 
is somewhat arbitrary as both mechanisms can occur si- 
multaneously without changing the coast fractality. Also 
a better model for damping should take care of a wave fre- 
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quency dependence as well as it should consider the exis- 
tence of localization by the frontier of the waves along the 
irregular coast 7]. This would modify Eq. (2) and change 
the time evolution. However, since percolation possess 
the universality properties of phase transitions 10], the 
fractal dimension of the coast should not depend on these 
factors. In the context of corrosion dynamics 15], this 
has been shown through arguments from dynamical field 
theory of absorbing states If large-scale modifica- 

tions or correlations exist in the lithology properties, then 
the resulting geometry would be more complex. But our 
main result, namely the existence of irregular coasts as 
a result of a self-stabilisation mechanism, would remain 
correct even though the geometry would be more com- 
plex than that of critical percolation. We believe that it 
is in those terms that can be interpreted the results of 
several detailed study of the self-similarity properties of 

sea-coasts nana. 

A more radical change in the geometry is expected if 
sediments transport is taken into account. In our simpli- 
fied model, the sediments are supposed to be transported 
offshore and disappear, while in the real erosion pro- 
cess these sediments are partially transported along the 
coast. However our approximation is justified by some 
examples, like the North Coast of Oregon (USA), a lead- 
ing hedge coast dominated by erosion with very limited 
coastal deposits [l7j . 

One should also mention that it has been found very 
recently that the GP power law s ap ply even when the 
front is too narrow to be fractal [20(. This extends con- 
siderably the range of application of GP to rocky coasts 
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which are irregular, but not to a fractal range. 

This work has presented a minimal model for the for- 
mation of fractal rocky coast morphology. This model 
bears on the reciprocal evolution of the erosion power 
and the topography of the coast submitted to that ero- 
sion: The more irregularly eroded the coast is, the weaker 
the average sea erosion power. Note that this seems to be 
an empirically known effect used to build efficient break- 
waters that are based on hierarchical accumulations of 
tetrapods piled over layers of smaller and smaller rocks, 
in close analogy with fractal geometry 21]. The retroac- 
tion leads to the spontaneous formation of a fractal sea- 
coast with a fractal dimension Df =4/3. The fractal ge- 
ometry plays here the role of a morphological attr actor: 
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when submitted to such type of erosion, forgetting its 
initial morphology. Note that our model suggests that, 
on the field, the islands which have resisted to an ero- 
sion under a force larger that the final force f(tf), are 
stronger that the coast itself. This could be verified on 
the historical data of known seacoasts and neighbouring 
islands evolutions. The model reproduces at least qual- 
itatively some of the features of real coasts using only 
simple ingredients: the randomness of the lithology and 
the decrease of the erosion power of the sea. It is worth 
to be noted that the use of simple geophysical ingredients 
leads to an evolution towards a self-organized fractality 
directly related to percolation theory. 
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